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Abstract. Motivated by the existence of unstable modes of compact stars that eventually grow 
large, we study the bulk viscosity of dense matter, taking into account non-linear effects arising 
in the large amplitude regime, where the deviation jUa of the chemical potentials from chemical 
equilibrium fulfills > T. We find that this supra-thermal bulk viscosity can provide a potential 
mechanism for saturating unstable modes in compact stars since the viscosity is strongly enhanced. 
Our study confirms previous results on strange quark matter and shows that the suprathermal 
enhancement is even stronger in the case of hadronic matter. We also comment on the competition of 
different weak channels and the presence of suprathermal effects in various color superconducting 
phases of dense quark matter. 
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INTRODUCTION 

Bulk viscosity is an important mechanism for the damping of oscillations of compact 
stars [1, 2, 3, 4, 5] and therefore is a useful measure for differentiating possible phases of 
dense matter in their interior. At sufficiently low viscosity and high rotation rate, r-modes 
are unstable and can cause rapid spin-down of the star via gravitational radiation [6] if 
their growth is stopped by some non-linear mechanism. Finding the relevant mechanism 
is important because it determines the amplitude at which the r-mode saturates, and 
hence the rate at which it spins down the star. At low amplitudes the bulk viscosity is 
amplitude-independent, but since the r-mode is unstable its amplitude grows, and unless 
stopped by other mechanisms [7, 8, 9, 10] will quickly enter the "supra-thermal" regime 
where the bulk viscosity grows with amplitude [11, 12, 13], and may become large 
enough to stop the growth of the mode. 

This proceeding is based on our recent paper [14], where we discuss the large am- 
plitude behavior of the bulk viscosity in detail. It expands the analytic approximation 
methods and extends the analysis to other phases of dense matter. 

GENERAL FORMALISM FOR THE BULK VISCOSITY 

The bulk viscosity of a given form of matter is defined by the response of the system 
to an oscillating compression and rarefaction. This corresponds to an oscillation in the 
densities of all exactly conserved quantities. In compact stars the relevant quantity is 
typically the baryon number density, 



n*(r,/) = n*(r) + 5n*(r,f) = n*(r) + An*(r) sin (co?) , 



(1) 



where n* is the equilibrium density and An* is the amplitude of the oscillation. We 
assume An* <C n*. The bulk viscosity is then given by [1 1] 



C~-/-\ (2) 
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Where {de/dt) diss is the averaged energy dissipation rate per volume in the fluid due to 
the oscillation. 

Due to the driving density fluctuation the quantity /1a = L; A*i — L/ A*/ i s driven out of 
equilibrium and its re-equilibration leads to the bulk viscosity. We will assume that the 
latter arises from beta-equilibration of fermionic species and that the net rate takes the 
general form 

r (-) = _fr^ A ^i + f Zj (0yj . (3) 

where ./V determines the highest power of /i A arising in the rate. The sum merely reflects 
that phase space for the weak processes is either provided by softening the Fermi surface 
due to finite temperature or by shifting the chemical potentials from their equilibrium 

values, so that K = 2N. In terms of dimensionless variables (p = cot and (<p) = 
the chemical fluctuation then fulfills the differential equation, 
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d cos (q>)-fs/[ l + £z^ 2j , (4) 



with the prefactors of the driving term and feedback term given by 

d ^ , fs *f. (5) 
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Where B and C are susceptibilities that can be obtained from the equation of state of 
matter. Note that the feedback term involves both linear and non-linear parts which are 
controlled by a single parameter / and that its particular form is determined by the 
constants %j which parametrize the particular weak rate. The viscosity is then given by 



C = -^co^AtJo Jo W)) 7 )Vcos((p)J(p. (6) 



where (<p;d,f) is the periodic solution to eq. (4) (for details see [14]). 

In the sub-thermal limit where /i A T, corresponding to srf <1, the non-linear terms 
in eq. (4) can be neglected, leading to the following equation for the sub-thermal bulk 
viscosity 

C 2 f T K = 2coBfT K 
^ co 2 + (BfT*) 2 W co 2 + (BfT K ) 2 ' { ' 

As can be seen from this equation, the sub-thermal bulk viscosity is amplitude- 
independent. As long as the combination of susceptibilities C 2 /B does not vary too 



quickly with temperature, the sub-thermal viscosity has a maximum 



At sufficiently low temperatures and high frequency one has / <?C 1 , and therefore the 
feedback terms in eq. (4) are negligible and the viscosity in this intermediate regime can 
be evaluated as 



C 2 TT« ( » (2./+1) lift /C An\ 2 ^ 
fe 0> 2 ^ + £i 2/(7 + 1)1 Vr 



(9) 



which can be combined with the subthermal result for an analytic solution in both 
regions ^ = C< + 6(T - T max ) £~. 

In the supra-thermal limit where /i A 3> 7\ corresponding to 3> 1, only the largest 
power of srf is relevant in eq. (4) and therefore the viscosity scales in this limit as 

<10) 

and decreases again at very large amplitudes. 

The general result for the bulk viscosity can be written in the following form 

t = ^ ax ^(dJ) = £^S(d,f) (11) 

where the dimensionless function that includes the non-trivial parameter dependence 
is given by 

2 r 27t 

J(dJ) = — / s/(r,d,f)cos(<p)d<p (12) 
ltd Jo 

The generic form of the function y(d,f) is shown in the left panel of Figure 1 It has 
a global maximum value of 1, reached in the sub-thermal limit and a line of local 
maxima along a parabola in the d-/-plane. Thus the maximum value (8) of the sub- 
thermal viscosity is also the maximum in the general case and depends only on the 
equation of state, the density and the frequency but is independent of the weak rate. The 
weak rate influences, however, at what temperatures and amplitudes the local maxima 
are reached. On the right panel of figure 1, the physically interesting regions of the 
parameters d and / for different forms of matter discussed below are compared to the 
regions of validity of the analytic approximations discussed above. Although the analytic 
approximations cover large parts of the parameter space, non-linear effects nevertheless 
require a numeric solution for extreme parameter values. 



BULK VISCOSITY FOR DIFFERENT TYPES OF DENSE 

MATTER 

In this section we apply our results given above to various models of quark matter and 
hadronic matter. We study strange quark matter [15] both within a simple quark gas 



FIGURE 1. Left panel: The function J' arising in the general solution eq. (11). The plot is given for the 
case of hadronic matter with modified Urea equilibration but the generic qualitative form is independent 
of the considered form of matter. The function has a global maximum of 1 reached asymptotically for 
d—>0,f= 1 and a line of slowly decreasing local maxima along a parabola in the d-f plane. The shading 
of the surface denotes the size of the amplitude srf so that dark shades of grey represent the supra-thermal 
regime. Eq. (5) relates d and / to underlying physical parameters such as temperature T and amplitude. 
An amplitude increase (keeping all other variables fixed) results in a linear increase in the variable d as 
shown by the dashed curves. An increase in temperature changes the viscosity along a line shown by the 
solid curves. Right panel: The relevant ranges of the parameters d and / in the general parametrization 
eq. (11) when the underlying parameters are varied in the physically interesting ranges T e [10 5 , 10 11 ] K, 
An/ne [l(T 5 ,l],Coe [1,1000] Hz, n e [1/4,10] no are denoted by the the dotted boundaries for hadronic 
matter with modified and direct Urea reactions and strange quark matter with non-leptonic reactions. 
The analytic sub-thermal approximation eq. (7) is valid well to the left of the black line, whereas the 
intermediate linear approximation eq. (9) is valid well below the piecewise curve formed by the gray 
linear segment on the left and the corresponding linear segment for the respective form of matter. 



model of free quarks in a "confining bag" as well as a phenomenological parameteri- 
zation of the quark matter equation of state proposed in [16] where a new parameter c 
has been introduced which incorporates some effects of strong interactions between the 
quarks. For hadronic matter we consider in addition to a free hadron gas the well-known 
model by Akmal, Pandharipande and Ravenhall [17] which relies on a potential model 
that reproduces scattering data at nuclear densities. As a low density extension of the 
APR data we use [18, 19]. In order to to apply our general results to other equations of 
state, away from chemical equilibrium we use the simple quadratic parameterization in 
terms of the symmetry energy employed in [20] . The resulting strong interaction param- 
eters describing the response of the different models are given in Table 1 . 

In unpaired strange quark matter, consisting of u, d and s quarks, the dominant channel 
for beta equilibration is the non-leptonic flavor changing process, d + u <-> s + u, where 
the equilibrating chemical potential in response to the the driving oscillation is given by 
Ma = M? ~ Hd- I n me case °f hadronic matter we assume that weak equilibration occurs 
via the Urea channel 



p + e — > n + v e , n — >■ p + e + v e 



(13) 



TABLE 1. Susceptibilities B and C for different model equations of state. In the 
case of hadronic matter with baryon density n a quadratic ansatz in the proton 
fraction x parameterized by the symmetry energy S is employed. The expressions 
for a free hadron gas are given to leading order in n/m^, and for quark matter with 
quark chemical potential \i q to next to leading order in m s /jj, g . The parameter c 
takes into account interaction effects within the employed quark matter model and 
vanishes for an ideal quark gas. 
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There are two qualitatively different cases depending on whether the direct process is 
possible or only the modified version where a bystander nucleon is necessary to satisfy 
energy-momentum conservation. Here the driving baryon number density oscillation 
yields the oscillating chemical potential difference \x^ = \x n — \x p — \x e . The equilibration 
rate parameters f , K and %u defined in eq. (3), for the weak processes discussed above, 
are given in Table 2. 



TABLE 2. Weak interaction parameters describing the considered damping pro- 
cess. Here fi q is the quark chemical potential, n is the baryon density, no nuclear 
saturation density and x the proton fraction. 
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Using eqs. (7-11) and the parameters given in Table 1 and 2, we obtain the results 
shown in Fig. 2, where the bulk viscosity as a function of the density oscillation ampli- 
tude at two temperatures is given. The solid lines show the results for interacting mat- 
ter whereas the dashed lines show the free hadron/quark gas results. The dotted lines, 
which are in most places invisible beneath the solid lines, show the combined analytic 
approximation below eq. (9) valid at moderate amplitudes. At the lower temperature 
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FIGURE 2. Comparison of the bulk viscosity of the different forms of matter studied in this work as a 
function of the density oscillation amplitude An/n. The frequency is co = 8.4 kHz, corresponding to an 
r-mode in a millisecond pulsar and n = 2«o. Left panel: low temperature T — 10 6 K; right panel: high 
temperature T = 10 9 K. The dashed curves are for the free hadron and free quark models; the solid curves 
are for APR hadron matter, and interacting quarks with m s = 150 MeV and c = 0.3 . 



the viscosity reaches the supra-thermal regime, already for small amplitudes, whereas 
at the higher temperature the sub-thermal regime, in which the viscosity is independent 
of the amplitude, extends to large amplitudes. The stronger non-linear feedback in the 
hadronic cases leads to a significantly steeper rise that correlates with the largest power 
in eq. (3). Interestingly, despite these differences the maximum value reached by varying 
the amplitude is still roughly the same as the maximum value in the sub-thermal limit 
eq. (8). This means that oscillations are approximately equally damped at all tempera- 
tures once the amplitude becomes sufficiently large. The maximum arises for amplitudes 
of the order 0.01, 0. 1 and 1 for strange quark matter and hadronic matter with direct and 
modified Urea, respectively. Therefore the simple analytic expression eq. (9) provides 
a remarkably good approximation for modified Urea processes. The supra-thermal en- 
hancement of the bulk viscosity is so strong that it could well provide the main saturation 
mechanism for unstable r-modes, stopping their growth at amplitudes that are below the 
threshold for other competing saturation mechanisms but large enough to allow spin- 
down of a neutron star via gravitational radiation on astrophysical time scales. 

The discussed suprathermal effects should be important for the competition of dif- 
ferent weak channels contributing to the viscosity. This has been considered for the 
non-leptonic and Urea channels in strange quark matter in [21]. From the stronger tem- 
perature dependence of the corresponding quark Urea rates it is clear that analogous to 
the direct hadronic Urea case the semi-leptonic contributions to the viscosity also have 
a stronger amplitude dependence than the non-leptonic contribution. Therefore, at large 
amplitudes the mixing between these two channels should be far more pronounced than 
what has been observed previously in the subthermal case [21]. Similar conclusions hold 
for the mixing in hadronic matter when the hyperonic channel is included [22]. 



Our general results for the bulk viscosity can also be applied to other forms of 
dense matter. For example in 2SC color- superconducting quark matter, where all strange 
quarks as well as up and down quarks of one color remain unpaired [23], it has been 
shown that at T <C T c the rate of the non-leptonic flavor changing process d + u -H- 
s + u differs from that in the unpaired quark matter only by a constant factor, T 2SC = 
T unpaired/ 9 [24]. The susceptibilities B and C receive additional contributions involving 
the gap parameter A, but since A < m s there is not a significant change and one can 
expect the same qualitative behavior for the bulk viscosity of the 2SC quark matter as in 
unpaired quark matter. According to eq. (11) the maximum value of the bulk viscosity 
Cmax is independent of the rate and thereby only slightly changed, but the bulk viscosity 
of the 2SC quark matter reaches its maximum value at higher amplitudes compared to 
unpaired quark matter. A similar statement holds for the bulk viscosity of ISC phases 
[25, 26]. Here the corresponding suppression factor depends on the particular pairing 
pattern, which may involve only part of the quarks on the Fermi surface, and can be as 
low as 0.03, thereby enhancing the effect discussed for 2SC pairing. 

In color-flavor locked (CFL) matter [23] all quark species are gapped and the rates for 
processes involving either quarks or massive pseudo Goldstone bosons, as well as the 
corresponding viscosities, are exponentially suppressed below the critical temperature of 
the order of tens of MeV [27]. The only ungapped modes are the superfluid phonons, but 
since their rate is not given by weak interactions [28] a simple estimate employing our 
general expressions is not possible in this case. Yet, we expect that neither the suppressed 
modes nor the phonons are very sensitive to changes in the underlying quark chemical 
potentials. Correspondingly, the amplitude dependence of the CFL viscosity should only 
be moderate and should not qualitatively change the picture obtained in [27, 28]. In 
contrast, in the kaon condensed CFLKO phase there are truly massless Goldstone bosons 
that mediate weak interactions [29]. Even though these particles are bosonic their weak 
reactions should be similarly influenced by changes in the corresponding flavor chemical 
potentials and therefore we expect a suprathermal enhancement at large amplitudes. 

CONCLUSION 

In this paper we have studied the bulk viscosity of dense matter taking into account non- 
linear effects that arise in the supra-thermal regime. We gave an approximate analytic 
solution as well as a general numerical solution for the bulk viscosity which are valid for 
all types of matter where equilibration occurs via fermions. We confirm previous results 
for the amplitude dependence of the bulk viscosity of strange quark matter [11] and 
find that these supra-thermal effects are parametrically even more important in nuclear 
matter. Whereas our result directly extend to partly gapped phases like ISC and 2SC, 
in the fully gapped CFL phase the suprathermal enhancement should be mild. The most 
obvious application of our results is to the damping of unstable r-mode oscillations in 
neutron stars. As the amplitude of the mode enters the supra-thermal regime the viscosity 
will increase steeply above the sub-thermal result and can exceed it by many orders of 
magnitude. It could therefore provide an adequate saturation mechanism if it occurs 
before other non-linear dynamic effects [7, 8, 9, 10] start to operate. 
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